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ABSTRACT 

Since fifteen years, the concentration-mass relation has been investigated diffusely in theoretical studies. On 
the other hand, only recently this relation has been derived from X-ray observations. When that happened, the 
results caused a certain level of concern: the X-ray normalizations and slopes were found significantly dissimilar 
from those predicted by theory. 

By analyzing a total of 52 objects, simulated each time with different physical recipes for the baryonic com- 
ponent, as well as 60 synthetic X-ray images, we aim at determining if these discrepancies are real or artificial. 
In particular, we investigate how the simulated concentration-mass relation depends z) on the radial range used to 
derive the concentration, ii) on the presence of baryons in the simulations, and on the prescription used to repro- 
duce the gas. Finally, we evaluate Hi) how the results differ when adopting an X-ray approach for the analysis and 
z'v) how the selection functions based on X-ray luminosity, temperature, and SZ-signal can impact the results. All 
effects studied go in the direction of explaining the discrepancy between observations and simulations, but they 
contribute at different levels: while the fitting radial range and the baryonic component play only a minor role, the 
X-ray approach and selection function have profound repercussion on the results. 

Extending the fit to the most central regions or reducing the fitting radius from the virial boundary to the typ- 
ical X-ray external radius causes an increase of the normalization. In the second case, a steeper slope is also 
found. Modifications in normalization or slope are more acute when baryons are included. Radiative simulations 
with only SN feedback present 30% higher concentrations than the DM case. The discrepancy is cancelled when 
the more powerful feedback by AGNs in considered. The concentration-mass relation derived from the X-ray 
synthetic catalogue is significantly steeper because the simultaneous influence of multiple effects such as environ- 
ment, dynamical state, dynamical history, mass bias, temperature bias, and their dependences on the radius and 
on the mass of the system. Finally, selecting the system via their SZ signal or mass-weighted temperature does 
not create any extra bias while a X-ray luminosity-based choice produces a net steepening of the relation since in 
the group-regime the most luminous systems are also the most concentrated. 

Subject headings: clusters: general — clusters: ICM — clusters: cosmology — X-rays: clusters 



1. INTRODUCTION 



The dark sector (dark energy and dark matter - DM) exceeds 
by about twenty-five times the baryonic content of our Uni- 
verse (e.g. |Voit||2005l |Borgani & Krav tsov 2011 Allen et al. 



2011 for recent reviews on cosmology with focus on galaxy 
clusters). As consequence, at first approximation, we may as- 
sume that simulations with only dark components properly de- 
scribe the salient halo properties such as the mass distribution, 
the tridimensional halo shape, the density profile, and the con- 
centration (see Kravtsov & Borgani|2012 for a galaxy-cluster 



formation review). However, an increasing number of theoret- 
ical papers begun to contest this simplistic view showing that 
baryons (in form of both hot gas and stars) can significantly 



influence all these topics (e.g. Rudd et al. 2008 


Stanek et al. 


2009 


Duffy et al.|2010 Cui et al. 2012[ Lau et al. 


2011 


Gnedin 


et al. 


20111 |van Daalen et al. 2011j|Zhu & Pan|12012; Mac- 



cio et al.|2012||Governato et al.|2 012; Bala guera-Antolinez &| 
Porc iani|2012| among the most recent publications on the var- 
ious aspects). Some disagreement between some observational 
results and theoretical predictions based on DM-only simula- 
tions support this discomfort. In this respect, a recent debate 



about the X-ray observed and simulated concentration-mass re- 
lation warmed up (Fedeli 2012 and references therein). Indeed, 
while earlier works based on 10-12 clusters observed by Chan- 
dra HVikhlinin et al.|2006]l and X MM-Newton ( |Pratt & Arnaud] 



2005 1 Pointecouteau et al.|2005 1 showed consistency with the 



oretical works (Dolag et al. 2004), lately, with the enlargement 
of the observational collections, more differences have arisen. 
The concentration-mass relation from the larger X-ray samples 
of |Buote et ap (120071 B 07), |Schmidt & Allenl ( |2007l SA07), 
and lEttori et al. d2010| E10) are significantly steeper (Fedeli 



|2012|l than the relation derived from the Millenium simulation 
(Gao et al. 2008, with results similar to aforementioned work 
by Dolag et al.|2004]>. A n exception to this mismatch is the re- 
sult by |Host & Hansen] ( |201 \\ who, however, analyzed only 11 
systems. 

The question we want to answer in this paper is whether 
the differences between theory and observations are real - i.e. 
some fundamental physical ingredients are missing in theoreti- 
cal models - or artificial - i.e. the comparison is conducted in 
an unfair fashion without considering possible biases that de- 
pend on the methods employed. 

All the mentioned papers on the c-M relations consistently 
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adopt the Navaro-Frenk- White (NFW) analytic model to fit the 
density profile and derive the concentration. Navarro et al. 
([1996} showed the existence of an universal profile that well de- 
scribes the halo density and mass distribution for a large range 
of halo mass and cosmologies ( |Navarro et al.|199"7] l. Their fit- 
ting formulae are characterized by two parameters: the normal- 
ization, p s , and the scale radius, r s . The expressions for the 
profiles of the NFW density and the NFW mass, as function of 
x = r/r s , are respectively: 



p{x)- 



(1) 



x(1+jc) 2 ' 

M«x) = 4wp s r 3 s f x , f x = Wl+x)-x/{l+x). (2) 
The logarithmic slope of the density profile measured at the 
scale radius is equal to -2, transitioning from the central asymp- 
totic value of-1 to the external value of -3. The normalization 
factor, p s , is connected to the characteristic overdensity, A, and 
to the critical density at the halo redshift z, p cr , through: 

Ap cr C 3 



(3) 



3 ln{l + c)-c/{\+c) 
where the critical density can be expressed as p cr = 3H 2 (z) / (8G), 
with H(z) the Hubble's parameter at redshift z and G the New- 
ton's gravitational constant. Finally, c is the concentration of 
the halo and it is defined as the number of times the scale ra- 
dius is contained within a fixed overdensity radius: c = R&/r s . 
In the following, we will always define the concentration in re- 
lation to the overdensity A = 20(Q 

The origins of the density shape mostly resides on the den- 
sity and triaxiality of the original peak and on the continuos 
contraction of the innermost material due to subsequent accre- 
tion during the collapse (Da lai et"aL"|2010| . This last point is 
also a suitable explanation for the increase of the concentra- 
tion parameter with the decrease of mass (e.g., NFW, Eke et al.] 



200T]|Bullock et al.|2001[|Dolag et al.|2004[ |Zhao et al.|2009| l 



The physical and cosmological justification is linked to the hi- 
erarchical scenario of structure formation: small halos form at 
earlier time when the universe was denser than today and, over 
time, the sedimentation of material at the center happens in an 
already-established high-density peak. 

The measure of a negative slope from observations is, there- 
fore, expected. The surprise is related to the values: the X-ray 
slope is up to 5 times steeper than what found in DM cosmo- 
logical simulations. In this paper, we want to test if this result is 
produced by hidden biases related to the two methods followed 
to analyze simulate and observed data. Indeed, it is unquestion- 
able that the basic procedural assumptions of the two analyses 
differ in several aspects. The simulated concentration-mass re- 
lation is mostly derived by fitting the NFW profile to the three- 
dimensional density distribution extracted considering the in- 
fluence of all the particles without exclusion of sub-structures 
or sub-clumps, and from the really central region to the out- 
skirts of the halo, usually up to the virial region, sometimes 
beyond. The sample selection is typically volume-limited with 
a well defined cut in mass. Finally, as said above, many theoret- 
ical works are based on large cosmological boxes of DM only 
(e.g[Bhattacha rya et al.|201 l||Kwan et al.|2012| among the most 
recent ones). On the contrary, the X-ray concentration-mass 
relation is derived from projected data that have a limited ra- 
dial range because of the signal-to-noise and the field of view. 

1 ^200 i s the radius of the sphere enclosing a mean density 200 times the critical 
density. Conversions of the NFW mass defined for d ifferent values of A are 
provided in the Appendix of Hu & Kravtsov ( 2003 1 and Appendix of Ettori 



Few clusters have been observed to 7?5oo with enough photon 
statistics and even fewer at /?200- For example, E10 (one of the 
samples we are comparing to) reaches on average 40% of /^oo 
(see Section |3.3.1| l. The observational selection function is not 
done in mass. In the best case, depends on the X-ray luminos- 
ity. In the worse (and the most common) situation, the clusters 
are chosen among the ones available on the archives, as conse- 
quence the selection function is not well defined. As final note, 
the real universe has baryons. 

The possibility that baryons might influence the structure- 
formation process and, in particular the DM distribution, has 
been introduced almost three decades ago with the develop- 
ment of the analytic model of 'adiabatic contraction' (|Barnes 



& White|1984||Blumenfhal et al.|1986||Ryden & Gunn|1987| i. 
More recently cosmological simulations push further the im- 
provement of the model (e.g. |Gnedin et al.]|2004| |G ustafsson 
et al. 120061 1 Abadi et al.|2"0T0l|Tissera et al.|2010||Gnedin et aL 
201 lj |Zemp et al.||2012| l. The underlying idea is that, at the 
center, the DM feels the growth of the total potential well due 
to the extra baryonic material accreted in consequence of cool- 
ing and condensation. However, even if a specific great effort 
has been devoted to study the concentration parameter in non 



Tisseraetal. 2010 



radiative simulations (e.g. Rasia et al. 2004; Pedrosa et al. 2009 



only recently these investigations have been 
extended to include various model of the intra-cluster medium 



(ICM) (|Rudd et al.||2008| |Duffy et al.pOlOl |Meneghetti et al. 
2010| |Macci6 et al.|20T2t|Governato et al.|2012[|Martizzi et al. 
20 1 2\ [Ragone-F igueroa et al. 2012! and also different dark- 
energy Trio^eTi~^Grossr^Springel 2009; B aldi|2012 De Boni 
etaL|20131 l. 

Summarizing, we aim at evaluating how the different steps 
and assumptions behind the two analyses can explain the dis- 
similar findings on the simulated and the X-ray-observed c-M 
relations. In particular we focus on the impact of z) the radial 
range, if) the presence of baryons, Hi) the X-ray analysis, and 
rv) the choice of the selection function. To reach these goals, we 
perform an intrinsic analysis on four sets of 52 clusters simu- 
lated each time with a different recipe for the ICM (runs include 
increasing level of complexity: DM only run, no-radiative ICM, 
cooling and feedback by mainly supernovae la and II, feedback 
by mainly active galactic nuclei). Finally, we follow an X-ray 
approach on the syn thetic Chandra-like catalogue introduced by 



Rasia et al. 



( |2012a|b| R12f] 



Let us restate and stress that the goal of the paper is not to 
provide another prediction for the c-M relations or to tackle 
the long-standing 'cusp-versus-core debate' (e.g. Ogiya & Mori 
Maccio et al.||2012| |Teyssier etlQT|2lIl3j |Martizzi etaL 



2011 



2012 i or the problem of the goodness of the NFW to describe 
density profiles of galaxy clusters (Meneghetti & Rasia, 2013). 
For these purposes, the size of our sample, even if comparable 
to observational works, is too small. We intend, instead, to 
check all procedural steps that can convey systematic biases 
affecting the c-M relation. We will discuss further the limits 
of our sample and possible influence on the results obtained. 

This paper is organized as follows: after presenting the sim- 
ulations in Section 2, we discuss the results derived from the 
intrinsic analysis of the simulated clusters in Section 3. The 
observational approach applied to synthetic X-ray images will 

A video abstract of R12 Rasia et al. |2012b| can be found at the following 
link: http://iopscience.iop.org/1367-2630/14/5/055()18?fromSearchPage=true 
The images of the synthetic catalogue can be found at http://www- 
personal.umich.edu/~rasia/ 



X-ray c—M relation 
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be presented in Section 4, while in Section 5 we will focus on 
the choice of the selection function. Our conclusions are out- 
lined in Section 6. 

2. SIMULATIONS 

In this Section, we briefly describe the four sets of simulated 
clusters that are used in our intrinsic analysis (Section 3). These 
four sets are obtained starting from the same initial conditions 
and using four different physical descriptions of the processes 
detemining the evolution of the baryonic component. More de- 
tails on the generation of the initial conditions are provided by 
|Bonafede efaLlpOlT) , while we refer to |Killedar eTaE] ( |2012l ) 
for a description of the different physical models implemented, 
and to |Rasia et al.| ( |2012a|b[ l for the characterization of the syn- 
thetic X-ray catalogue of X-ray. 

2.1. The sets of simulated clusters 

Our simulated clusters have been chosen within twenty nine 
Lagrangian regions identified within a 1 h~ 3 Gpc 3 low-resolution 
N-body cosmological simulation^] The cosmological model as- 
sumed is a flat ACDM model, with Q. m = 0.24 for the matter 
density parameter, Obar = 0.04 for the contribution of baryons, 
Ho = 72kms" 1 Mpc" 1 for the present-day Hubble constant, n s = 
0.96 for the primordial spectral index and erg = 0.8 for the nor- 
malization of the power spectrum. Resolution is increased within 
the selected regions by enlarging the number of particles and 
correspondingly adding higher-frequency modes from the power 
spectrum of the same cosmological model, by using the Zoomed 
Initial Condition (ZIC) technique (Torme rTet al.|[1997| . The 
runs were carried out using the GADGET-3 code, a newer and 



more efficient version of the former GADGET-2 code ( Springel 



2005 ). In all simulations a Plummer-equivalent softening length 



for the computation of the gravitational force in the high-resolution 
region was fixed to e = 5 /r'kpc in physical units for most recent 
redshift (z < 2), while was kept fixed in comoving units at ealier 
epochs. As for the computation of hydrodynamic forces, we as- 
sume the SPH smoothing length to reach a minimum allowed 
value of 0.5e. 

Besides a set of dark matter only (DM) simulations and a 
set of non-radiative (NR) hydrodynamical simulations, we also 
carried out two sets of radiative simulations (CSF and AGN ). 
In these simulations radiative cooling rates are computed by 
following the procedure presented by Wiersma et al. (2009), 
including the effect of CMB radiation and of UV/X-ray back- 
ground radiation from quasars and galaxies ( |Haardt & Madau 
2001 ). The contributions to cooling from eleven elements (H, 
He, C, N, O, Ne, Mg, Si, S, Ca, Fe) have been pre-computed 
using CLOUDY ( |Ferland et al.|1998] ) for an optically thin gas in 
(photo-)ionisation equilibrium. 

The CSF set of simulations includes star formation and the 
effect of feedback triggered by supernova (SN) explosions. As 
for the star formation model, gas particles above a given thresh- 
old density are treated as multiphase, so as to provide a sub- 
resolution description of the interstellar medium, according to 
the model originally described by Springel & Hernquist (2003 ). 
SN-II, SN-Ia and low and intermediate mass stars contribute to 
the production of metals according to the model described by 
Torn atore et al.| ( |2007[ ): assuming a Chabrier IMF (Chabrier 
2003 ) stars produce metals over the time-scale determined by 

" 1993] ). 



The effect of kinetic feedback triggered by SN-II is included ac- 
cording to the prescription by Springel & Hernquist ( 2003| ). In 
the CSF simulation set we assume v w = 500km s _1 for the ve- 
locity of the outflows, with a mass-upload rate that is two times 
the value of the local star formation rate. 

The AGN set of simulations is carried out by including the 
same physical processes as in the CSF case, with a lower wind 
velocity of v w = 350 km s" 1 , but also including the effect of 
AGN feedback. In this model, largely based on the original im- 
plementation of black-hole (BH) feedback by Spri ngel et al.| 
(2005), energy results from gas accretion onto supermassive 
black holes. BHs are included as sink particles, which grow 
in mass by gas accretion and merging with other BHs. Gas 
accretion proceeds at a Bondi rate, and is limited by the Ed- 
dington rate. Once the accretion rate is computed for each BH 
particle, a stochastic criterion is used to select the surrounding 
gas particles to be accreted. Unlike in Springe l et al.|2005] in 
which a selected gas particle contributes to accretion with all 
its mass, we included the possibility for a gas particle to accrete 
only with a slice of its mass, which corresponds to 1/4 of its 
original mass, thus providing a more continuous description of 
the accretion process (see also Fabja n et al.||2010] ). BH parti- 
cles are initially seeded with a mass of 5 x 10 5 h _1 MQ. Seeding 
of BH particles takes place at the minimum of the potential of 
halos when they first reach a minimum friend-of-friend mass of 
2.5 x lO 13 h _1 M0 (using a linking length of 0.16 in units of the 
mean interparticle separation), with the further condition that 
such halos should contain a minimum mass fraction in stars of 
0.02. This condition guarantees that substantial star formation 
took place in such halos and that seeding does not take place in 
structures linked together by the friend-of-friend algorithm, not 
corresponding to a real halo. A radiation efficiency parameter 
of e r = 0.1 determines the amount of radiated energy extracted 
from gas accretion, in units of the rest-mass energy of the ac- 
creted gas, with the BH mass being correspondingly decreased 
by this amount. A fraction ej of this radiated energy is ther- 
mally coupled to the surrounding gas. We use ey = 0.05 for this 
feedback efficiency, which increases to e/ = 0.2 when accretion 
enters in the quiescent "radio" mode and takes place at a rate 
smaller than one-hundredth of the Eddington limit (e.g. Sijacki 
et al.|2007||Fabjan et al.|20T0l ). 



the mass-dependent life-times of Padovani & Matteucci 

3 h = 0.72 throughout the paper 



The analysis described in Section 3 is based on a sample of 
52 objects. These are the central clusters of each of the 29 re- 
gion plus three halos selected from rich environment (regions 
containing more than 20 halos with virial mass greater than 
5 x IO^/t'Mq), six from poor environment (less than 8 ha- 
los per region in the same mass limit), and fourteen lying in a 
medium populated environment. The covered mass interval is 
similar to that usually considered in observational works, with 
M 2 rjo ranging from 7.5 x lO 13 ^ 1 M to 2.5 x lO^h" 1 M . 
All clusters have been analyzed at redshift zero (for comparison 
with other works present in literature) and 0.25 (for comparison 
with the X-ray sample of R12). Throughout the paper, we will 
focus on the z = objects reporting our z = 0.25 results in the 
Appendix. 

2.2. Synthetic catalogue 

The X-ray synthetic catalogue employed in Section 5 and 6 
contains 60 event files related to three orthogonal projections of 
twenty massive objects part of the CSF sample at z = 0.25. 

The X-ray catalogue has been created using the X-ray MAp 
Simulator (X-M AS |Gardini et al.|2004]|Rasia et aLp0 08) after 
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the removal of the over-cooled particles identified in the den- 
sity (^[gr/cm 3 ]) - temperature (r[keV]) plane that satisfy to 
the condition: T < 3 x 10 6 p° 25 . This cut reduces the presence 
of small peaked clumps in the X-ray images without affect- 
ing the overall distributions of density, temperature, and X-ray 
luminosity (see Appendix of R12). The redistribution matrix 
function and ancillary response function adopted are those of 
Chandra ACIS-S3. The field of view is 16 arcmin, equivalent 
to ~ 2.5 Mpc for our cosmology and redshift, and the exposure 
time is of 100 ksec. The metallicity is kept constant to 0.3 solar 
with respect to the tables of Anders & Grevesse ( 1989[ l with a 
small correction on the helium abundances to be consistent with 
the simulated hydrogen mass fraction. The galactic absorption, 
described by a WABS model, is fixed to Nh = 5 x 10" 20 cm -2 . 
Along the line of sight, we consider the information of all the 
particles located up to 5 h _I Mpc (in both directions) from the 
center of the object. The final event files avoid the inclusion 
of the background since this component does not influence the 
mass profile derivation (Rasi a et al.|2006 2012b I. 



SIM radial range 



3. INTRINSIC C-M 

3.1. Density profiles and fitting procedure 

To perform the tri-dimensional analysis of our simulated sam- 
ple, we proceed by extracting the spherical mass profiles in 50 
bins logarithmically spaced from 10 h" 1 kpc to 5 h" 1 Mpc and 
centered on the miminum of the potential well. For each clus- 
ter, all the bins external to s» 2% of Z?200 have more than one 
thousand particles. This number represents the threshold for 
numerical convergency of the inner slope ( Power et al.||2003' 



Following the observational approach, we fit our simulated pro- 
file with the NFW mass expression (Eq. 2) to recover the two 
parameters p s and r s . The process is performed multiple times 
on the same object changing in each circumstance the radial 
range assumed (Section [3.3.1[ l. As subsequent step, we derive 
the NFW concentration using Eq. [3 The concentration-mass 
relation is, then, calculated by linearly fitting the following ex- 
pression: 

i i i M NFW 

log 1()C = log 10 c„ + a x log 10 5xlQl4h _ lM0 (4) 

For the hydro-simulations (NR , CSF , and AGN ), we analyze 
both the total mass profile and the profile of the only dark mat- 
ter. In each radial bin, the NFW fitting procedure assumes an 
error on the mass equivalent to 10% of the mass itself, value 
consistent with observations and previous analyses on X-ray 
synthetic catalogues (Meneghetti et al. 2010, R12). 

3.2. Redefinition of the sample 

When studying the influence of the radial range (Section 3.3) 
and the baryonic physics (Section 3.4 1, we consider only the 



profiles that show a good consistency with the NFW descrip- 
tion. In this way we exclude both objects that presenting one 
or more large substructures produce a secondary peak in the 
density profile and halos that do not have a coherent concentra- 
tion parameter associated. This a pr/on'-selection is not usually 
embraced in theoretical works, however, it is often implicitly 
adopted in observational samples especially if they are focused 
on regular systems. The careful choice of the simulated clus- 
ters that need to be discarded is done on the basis of the NFW 
residuals: 

^n„. ™n i„„ e-KM vi2 

(5) 
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FIG. 1. — Distribution of residuals as defined in Eq|5]of the four physics 
and at redshift zero. The concentration is derive d assu ming the SIM radial 
range, [0.06- 1.27] X i?200> as defined in Section |3.3.1| Acceptable clusters 
have residuals below 0.035 (red dotted vertical line). 



For all situations considered (varying ICM physics or radial 
ranges) and for both redshifts, we study the distributions of the 
residuals and notice that the majority of the objects has resid- 
uals below 0.035 at z = (see Figure [T]) and 0.05 at z = 0.25. 
The residuals are significantly larger at z = 0.25 especially for 
the NR physics. Furthermore, at fixed redshift, the outliers are 
often the same objects despite the change of the simulated set 
or the fitting radial range. Per each redshift, we identify these 
systems: nine clusters at z = and eleven at z = 0.25. The NR 
physics presents larger deviations because accreted substruc- 
tures requires more time to virialize and thermalize under the 
no-radiative physical condition (Dola g"et al.|2009[ ). 

In the following, we decide to always exclude these identi- 
fied systems and any other existing outliers (defined as having 
a > 0.035 at z = and a > 0.05 at z = 0.25). For each situa- 
tion investigated, before deriving the c-M relation, we check 
whether the excluded halos are spread across the mass range. 
When this does not happen (for example, if all object below a 
certain mass disappear) we prefer to not show the derived c-M 
relation because of the significant change on the statistical con- 
ditions. 

3.3. Radial range 

3.3.1. Typical theoretical and observational radial range 

Previous works demonstrated that the fitting radial range has 
some level of impact on the derived NFW concentration. For 
example, Rasia et al. (2006) showed that the concentration of 
a single clusters varies from ~ 5.3 to 7 (23% variation) when 
the outer radius moves from the virial radius to 20-30% of that 



2 E [log io (Md - \og 10 (M NFW j)Y 
a = 



N hil 



radius. Using their X-ray observations on A2717, Gastaldello 
|et al.| ( |2007| l warned on the derivation of the concentration pa- 
rameter if the data do not extend beyond the measured scale 
radius. Fedeli ( 2012| l stressed more on the consequences of the 
choice of the inner radius: fitting from 1% of R 200 induces an 
higher concentration, especially in small systems, than starting 
from 5% of 7?2oo- 

To evaluate the influence of the radial range assumed to fit 
the mass profile with the NFW formula, we begin by studying 
the DM sample. We first consider radial ranges that are 'typi- 
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cally' used in theory and X-ray studies, then we move to a more 
general discussion. 

Our reference for theoretical works are Net o et aT] ( |2007| l, 
puffy et al.| ( [2008][20T0l > and Meneghetti & Rasia (2013). Their 
radial range used to fit the NFW formula is between 5% and 
100% of the virial radius. The innermost limit was set to satisfy 
the requirement on the numerical convergence (Section 3.1 and 
Power et al. 12003 ). As done by most observational works, we 
prefer to report all the measurements to /?200> that is equal to 
R V ir/\.21 for our sample, redshifts, and cosmology. Therefore, 
the typical simulation radial range is [0.06,1.27] x /?200- We 
label this radial range as SIM. 

The X-ray radial range is, instead, linked to the one used in 
Ettori et al. (2010 1. The Authors provided the radial bound- 
aries used to compute the NFW best fit. The inner radius was 
constantly fixed to 50 kpc to exclude the impact of the central 
galaxy on the density profile and thus avoiding the stellar com- 
ponent influence. In our whole sample, this value is larger than 
0.03 x /?200 fulfilling the convergency requirement. E10 also 
furnished two outer radii, one associated to the surface bright- 
ness profile and the other to the spectroscopic temperature pro- 
file. Per each cluster, we select the minimum of the two and 
compare its values to the derived /?5oo- The mean ratio resulted 
to be equal to 0.6^] Concluding, we label the radial range from 
50 kpc to 0.6 x /? 5 (x|^]as X-RAY. 

Considering both the SIM and the X-RAY radial ranges we fit 
the same DM cluster mass profiles, derive values of both NFW 
concentration and mass, and fit the concentration-mass relation 
with Eq. 4. For the z=0 sample the resulting normalization, Cq, 
and slope, a, are reported in the first two rows of the first sec- 
tion of Table [T] with their 1 a error. The two slopes are consis- 
tent within 1 (7 while the normalization, when a typical X-RAY 
radial range is used, increases by 6-10% with respect to the 
theoretical radial range. Limiting the fit to inner parts of the 
mass profile produces an increase of the halo concentration at 
all mass scale. The cause is the combined effects of lowering 
both the internal and external boundaries. In this circumstance, 
the innermost part of the profile weights very strongly. Indeed, 
within 100 h" 1 kpc almost all the NFW X-RAY best fits are 
lying above the SIM ones. Comparing the concentrations, the 
scale radii, and the two 1?200 derived, we found that with respect 
to the SIM radial range: i) the X-RAY Z?20o(= r s x c) is only 4% 
lower; ii) the X-RAY scale radii are significantly lower (~ 13%) 
and, consequently, Hi) the X-RAY concentrations are higher (~ 
16%). 

3.3.2. Radial range generalization 

To broaden the discussion about the radial ranges, we ex- 
plore some variation on the external and the internal radii with 
respect to the SIM radial range. The results are shown in Fig- 
ure [2] and listed in Table [T] For clarity, we show the cluster 
points only for the SIM radial range. In the top panel, the in- 
ner radius is constantly equal to 0.06 x/?200 while the external 
one varies from 0.4 x7?2oo (solid blue line) to 2 x R 200 (dashed 
orange line). From this plot, it appears that both slopes and nor- 
malization slightly increase when the external radius becomes 
smaller and closer to the X-RAY limit (0.6 x 7? 5 oo ~ 0.4 x 7?20o)- 

Even if the changes are very small and within 1 a, we tested 



4 In Ettori et al. 20101 two different R500 radii were reported accordingly to 
two different derivations of the mass profiles. For our purpose, the variation 
between these are minimal, as we tested. 

5 this value corresponds roughly to 0.4 X R200 f° r our objects. 



Table 1 

Concentration-mass best-fit parameters with 1 a error 

of the dm sample at redshift z = 0. Apart the 
x-rayradial range fixed between 50 kpc and 0.6 x /?500, 
all remaining radial ranges are defined in units of 
#200, WITH SIM=[0.06-1.27] xR 2 oo- 



fitting 
radial range 


DM 

N a co 


SIM 
X-RAY 


43 -0.156 ± 0.033 3.724 ± 0.015 

42 -0.114 ±0.034 4.014 ±0.015 


[0.06-0.4] 
[0.06-2.0] 
[0.03-1.27] 
[0.2-1.27] 


43 -0.153 ±0.034 3.844± 0.016 
43 -0.137 ±0.033 3.722 ± 0.015 

42 -0.123 ±0.032 3.827 ± 0.015 

43 -0.142 ±0.037 3.597± 0.017 


DM-only 
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FIG. 2. — Each line represent a different concentration mass-relation ex- 
tracted by varying the limits of the SIM radial range. In the top panel, we vary 
the external radius from 0.4 (solid blue line) to 2 X R200 (dotted-dashed orange 
line). In the bottom panel, the internal radius moves from 0.03 (solid red line) 
to 0.2x^200 (dotted-dashed green line). The points are the values of NFW 
concentrations and masses computed assuming the SIM radial range. The big 
crosses indicate the 9 clusters excluded for their large residuals. 



the origin of this trend. We found that halos with different 
masses react differently to changes of the radial fitting: smaller 
halos play the major role. For example, the parameters p s and 
r s , of our most massive clusters (M tlue > 10 15 h _1 M Q ^]do not 
depend on the radial range chosen implying an averaged small 
variation on the NFW masses and concentrations (< 5%). The 
largest contributors of the slope variation are objects with low 
mass. Our smallest clusters (M tlue < 1.5 x 10 14 h _I M ) show a 
change in the scaling radius and in the normalization of ~ 20% 
and ~ 25% when moving from the widest ([0.06-2] x R200) 
to the narrowest ([0.06-0.4] x R200) radial range. The con- 
sequence is to produce a ~ 10% larger concentration and 20% 
smaller NFW mass with a shift not parallel to the SlMconcentration- 
mass relation but almost orthogonal to it causing the overall 
change of the slope. We conclude that the radial range might 
play a more relevant role whenever the sample extends to small 

6 We indicate with M uue the true mass of the simulated objects within i?2oo- 
Notice that it might differ from the NFW masses reported in the various Figures. 
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halos. 

The other scenario presented in the lower panel of the same 
Figure[2]regards the change on the inner radius from 0.03 xR 200 
(dotted-dashed red line) to 0.2 x Z?200 (solid green line) while 
the external radius is kept fixed at 1.27x/?200- Excluding a sig- 
nificant portion of the inner profile (with a cut similar to that 
applied in weak-lensing studies) causes a reduction on the nor- 
malization. Vice versa, extending the fit to inner regions (i.e. to 
values comparable to the X-RAYinner radius) increases the nor- 
malizations as previously noticed comparing the X-RAY and 
SIM radial ranges. Also in this case, however, the modifications 
in the concentration-mass relation are not enough to explain the 
differences between the simulated and observed relations. 

3.4. Impact of baryons: c-M relation 

To evaluate the influence of baryons and of the baryonic physics 
on the c-M relation, we fit the NFW formula (Eq. 2) to the 
mass profiles of the hydro simulations using all the radial ranges 
discussed above (SIM, X-RAY, and the four cases of Figure [2] 
and Table 1). On top of the total mass profile, in this Section, 
we have considered also the DM-only mass profile, i.e. the 
mass profiles derived from the DM particles within the hydro- 
simulations (thus, ignoring both baryonic components: gas and 
star). 



SIM radial range 
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FIG. 3. — Profiles of the residuals for all the hydro simulation (NR in 
green, CSF in red, and AGN in magenta) at redshift zero. The NFW fitting 
is deriv ed ass uming the SIM radial range ([0.06- 1.27] X R200) as defined in 
Section [3.3.1) . In the top panel, the NFW fit is applied to the total mass profile, 
while in the bottom one to the mass profiles of the DM-particles. 



3.4.1. Deviation from the NFW formula 

The number of clusters whose mass profiles are well fitted 
by the NFW formalism (see Section |3.1| l is also listed in the 
same Table. The CSF sample is heavily reduced due to large 
deviation from the NFW formula especially at small radii. This 
is evident looking at the top panel of Figure [3] where we show 
the average profiles of the NFW-mass deviations from the true 
mass profile: (MNFw _ Mrae)/^true- The Figure built excluding 
the nine systems with large residuals (see Section [3~T) i. Even 
so, the CSF clusters (in red) largely depart from a pure NFW at 
the center. The NR clusters (in green) are very well represented 
by the NFW formula showing almost zero deviations at all radii. 
The AGN runs are halfway. Since they include radiative cooling, 



as the CSF runs, they also have central profiles steeper than an 
NFW profile and than the NR clusters. However, the presence 
of AGNs, which push a considerable amount of gas towards 
more external radii, reduces this effect with respect to the CSF 
simulations (Fabja net al.|2010) . 

The condensation of baryons in the center is very efficient in 
less massive objects. Indeed, almost all our CSF groups with 
M trae < 2 x 10 14 h _I M Q have residuals a > 0.035 when the fit is 
extended to a small innermost limits such as either 0.03 x R20Q 
or 50 kpc. The exclusion of all these systems implies a drastic 
reduction of the mass range impeding to robustly derive the c— 
M relation. 

From the bottom panel of Figure [3] we learn that the NFW 
is an excellent description of the mass profiles of the DM-only 
particles. Looking at similar profiles generated by varying the 
radial range adopted (not shown), we notice that this feature 
persists at all times. As further confirmation, when we restrict 
our analysis to the dark component only few objects show large 
residuals. Indeed, in Table [2] the values of N associated to the 
DM mass profiles within the hydro-simulations are almost al- 
ways equal to 43. 

3.4.2. Comparison with the DM set: slope 

The results of the baryonic simulations, presented in Table|2] 
indicate that the slopes of the c-M relations are always consis- 
tent within 1 a error with the DM slopes (Table 1). However, 
when baryons are included we see an amplification of the ef- 
fects discussed in Section 3.3. For example, enlarging the ra- 
dial range from /?,,,> (as in SIM) to twice /?200 has the effect 
of always reducing the c-M slopes by at least 30% (twice as 
much as for the DM runs). In this circumstance, filaments or 
substructures are entering in the scene and reduce the concen- 
tration of small systems. When baryons are present the influ- 
ence of the environment in more prominent for their collisional 
nature. The deviations from the NFW formula at large radii, 
in fact, are common to all three physics and are in agreement 
with the finding of [Dolag et al.| ([2006) who showed that halos 
in a filamentary structure or in the vicinity of massive systems 
have an external slope that differ from an NFW and it asymp- 
totically tends to -2 instead of to -3. In addition, reducing the 
external radius from 2 x /?20o to 0.4 x /?2oo produces a more sig- 
nificant steeping of the c-M relation: the two a values change 
by ~12% for DM runs, while the difference becomes ^40-45% 
for NR and CSF ~65% for AGN . 

Looking with more details at the AGN case, in Figure |4j we 
plot the averaged density profiles of the least massive (black 
line) and most massive (red line) clusters normalized at Rzca- 
The profiles are not exactly self similar: least massive groups 
are on average more concentrated than largest clusters. In the 
bottom panel, we show how much concentrations and NFW 
masses vary in relation to the external radius. Most points are 
lying in a precise location in this plot because the concentra- 
tions and masses of the two cases refer to the same potential 
well as we are modifying only the radial range of the fitting 
over the same mass profile. A certain degeneracy between the 
parameters is present. We find that our most massive systems 
(M true > 10 15 h~ l M Q ) behave exactly as in the DM simulations: 
their concentrations and masses vary very little with the change 
of the fitting radial range (the red squares in the central box 
show small deviations, less than 2-3%). On the other hand, 
the smallest systems (M tme < 1.5 x 10 14 /i _1 M Q ) show a drop in 
concentration and an increase in mass when we augment the 
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Table 2 

Concentration-mass best-fit parameters with 1 a error. The fitting radial ranges are defined in units of R 2 m (with 
siM=[0.06-1.27]i?2oo)- The only exception is the X-ray radial range defined between 50 kpc and 0.6 x R 5 qo- The block on 
the left corresponds to results obtained by fitting the total density profile while the block on the right refers to 

the results related to the dm only density profile. 







NR 






DM in NR 




N 


a 


co 


N 


a 


co 


SIM 


43 


-0.158 ±0.038 


3.978 ± 0.018 


43 


-0.165 ±0.037 


3.688 ±0.017 


X-RAY 


42 


-0.159 ±0.041 


4.000 ± 0.019 


42 


-0.153 ±0.038 


3.826 ±0.018 


[0.06-0.4] 


43 


-0.174 ±0.041 


4.045 ± 0.019 


41 


-0.172 ±0.039 


3.777 ±0.018 


[0.06-2.0] 


41 


-0.119 ±0.035 


4.079 ± 0.016 


42 


-0.136 ± 0.036 


3.754 ±0.017 


[0.03-1.27] 


41 


-0.132 ±0.038 


4.016 ± 0.017 


41 


-0.141 ± 0.035 


3.791± 0.017 


[0.21-1.27] 


43 


-0.139 ±0.039 


3.891 ±0.018 


43 


-0.158 ±0.039 


3.542 ±0.018 






CSF 






DM in CSF 




N 


a 




N 


a 


co 


SIM 


40 


-0.205 ± 0.038 


4.781 ±0.017 


43 


-0.166 ±0.035 


3.925 ±0.016 


X-RAY 








42 


-0.129 ± 0.033 


4.612 ±0.016 


[0.06-0.4] 


43 


-0.224 ± 0.034 


5.457 ± 0.016 


43 


-0.167 ±0.036 


4.194 ±0.017 


[0.06-2.0] 


36 


-0.159 ±0.038 


4.739 ± 0.016 


43 


-0.137 ±0.032 


3.968 ±0.015 


[0.03-1.27] 








39 


-0.190 ±0.031 


4.287 ±0.014 


[0.21-1.27] 


43 


-0.154 ±0.038 


4.042 ± 0.017 


43 


-0.160 ±0.038 


3.568 ±0.018 






AGN 






DM in AGN 




N 


a 


co 


N 


a 


co 


SIM 


43 


-0.130 ±0.037 


3.856 ±0.017 


42 


-0.155 ± 0.037 


3.631 ±0.017 


X-RAY 


37 


-0.078 ±0.042 


4.891 ±0.020 


43 


-0.107 ± 0.039 


3.925 ±0.018 


[0.06-0.4] 


43 


-0.149 ±0.04 


4.072 ± 0.019 


43 


-0.134 ±0.040 


3.748 ±0.019 


[0.06-2.0] 


41 


-0.090 ±0.034 


3.943 ±0.015 


41 


-0.098 ±0.034 


3.701 ±0.016 


[0.03-1.27] 


24 


-0.084 ±0.041 


4.113 ±0.019 


40 


-0.147 ±0.036 


3.745 ± 0.016 


[0.21-1.27] 


43 


-0.107 ±0.038 


3.641 ±0.018 


43 


-0.146 ±0.038 


3.406 ±0.018 
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outer radius (the black square in the central box are all consis- 
tently different from zero). 

3.4.3. Comparison with the DM set: normalization 

More generally, when we compare the c-M normalizations 
of the hydro-simulations with the DM case, we see a consis- 
tent increase (Table|2|. The CSF runs are those with the biggest 
deviation, especially when we reduce the external limit to 0.4 
Xi?200- This result is visualized in the top panel of Figure [5] 
where similarly to what shown in the lower box of Figure [4] we 
are reporting the variation in NFW concentrations and masses 
moving from the DM runs to the hydro-simulations using the 
SIM radial range and for z = 0. The CSF simulations produced 
a steeper relation but also an higher normalization. Indeed, the 
inclusion of baryons without a strong feedback induces an in- 
crease of the concentration at all mass scale (see lower panel in 
the Figure). The effect is mitigated, if not completely canceled 
when we include the strong feedback produced by the AGNs. 
From the bottom panel, we can evince that the slope of AGN 
and NR is very similar to the DM one (both sets of green and 
magenta points are parallel to the DM black horizontal line). 
Finally, we observe that the deviations in both concentration 
and mass are correlated (top panel) implying that even if the 
simulated underlying physics is changed the potential wells of 
galaxy clusters are not amply modified. 

AGN 
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FIG. 4. — Top panel: averaged mass profiles multiplied by the radius squared 
of the least (black line) and most (red lines) massive clusters of the AGN sam- 
ple. Bottom panel: differences between the NFW concentrations and the NFW 
masses: A Q4 = (c2-c4)/c2 and A M , 4 = (M2-M4)/M2, where c2 (c4) and 
Ml (M4) refer to values derived in the [0.06-2] X R 2 00 ([0.06-0.4] X i? 2 oo) 
radial range. 



3.4.4. Comparison with previous theoretical works 
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FIG. 5. — Top panel: variation in the NFW concentrations and in the masses 
when the quantities are computed in the hydro-simulations with respect to the 
DM runs. Red, green, and magenta triangles represent CSF, OVISC, and AGN 
results from the fitting in the SIM radial range, respectively. Bottom panel the 
only deviation of the concentration in function of the cluster mass. 



Rudd et al. ( 2008| ), studying the influence of baryons on the 
matter power spectrum, analyzed how baryons affect density 
profiles and concentrations. Their results are very similar to 
ours: non radiative simulation have a concentration boosted 
by 5-10% while simulations with gas cooling, star formation, 
and SNe feedback present a larger difference (between 20 and 
50%). Similar values have been found by De Boni et al. ( 201 3) ) 
with some variations caused by the particular dark energy model 
considered. Duffy et al. ( 2010| l studied, as we do, also runs 
with the stronger feedback by AGN. Unfortunately, they re- 
ferred their studies to smaller masses, thus, a straightforward 
comparison cannot be done on the whole sample. However, 
their most-massive bin (M ~ 10 14 h" 1 M ) returns a concentra- 
tion very close to their DM simulation as it is in our case. 

3.4.5. Comparison with previous observational works 

To compare with B07, SA07, and E10 we will use the re- 
analysis performed by FedeTI|(2012i who to homogenize the 
results, i) referred to an overdensity of 200 and to the latest 
WMAP results Komatsu et al. (2011) ; ii) removed the relation 
redshift dependence, Hi) chose a fixed pivot point at mrj = 5 x 
lO^h-'M©, as we do. 

We refer to all these papers for a complete description of their 
samples, here, we summarize the salient points: 
B07 studied 39 objects with masses from 6 x 10 12 M Q to 2 x 
10 15 M Q . The core of their sample resides in early-type galax- 
ies and groups (16 systems of their sample have mass below 
10 14 M Q ) to which they added the massive clusters from the 
sample of Pointecouteau et al. ( 2005} and |Vikhlinin et al.| (f2"006 ). 
Restricting their analysis to only clusters with mass above lD~ pr M© 
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B07 confirmed what previously found by these latter works. 
SA07 used 34 massive halos having mass-weighted temper- 
atures within /?2500 above 5 keV and showing regular X-ray 
isoflux contours and minimal isophote centroid variations. Most 
of the X-ray information are contained in /?2500 (~ 0.3 x /?20o)- 
Finally, E10 considered 44 very luminous nearby clusters with 
a minimum mass of 1O 14 M0. 

In Figure[6] the CSF simulations are shown in red, NR runs in 
green, and AGN ones in magenta. For each physics, we chose to 
show the values obtained from the SIM radial ranges with thick 
solid lines in the case of the total density and with thin solid 
lines for the DM component in each physics. The vertical and 
horizontal lines and the black region at their intersections repre- 
sent the values of the DM simulations (bold values of Table [TJ. 
Finally, the larger dotted-dashed rectangular regions indicate 
the location in the parameter space identified by varying the ra- 
dial ranges (i.e. values between the minimum and maximum 
listed in Table [2]). 

The figure clearly confirms our findings described in the pre- 
viously and summarized in the Table [2j variation in the gas 
physics do not lead to big changes in the slope. The highest 
values, registered for CSF clusters, is 40% higher than the DM 
slope. However, it is still much smaller than the values found 
by |Schmidt & Allen|p007l l and |Ettori et alJpOlO i. 

The normalization, instead, is affected by a larger degree. 
We notice that for this parameter not only the CSF clusters but 
also those with AGN feedback are in agreement with |Buote| 
|et al.| ( |2007| l when we consider the X-ray radial range. An- 
other agreement is between the normalization of the DM com- 
ponent and the value derived by Ett ori et al.| ( |2010| > who also 
analyzed the DM component. Our agreement with E10 persists 
also when we compare with their normalization derived from 
the whole sample and after fixing the slope equal to a = -0. 1 . 
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FIG. 6. — Normalization, cq, and slope, a, of the c — M relations derived in 
the X-ray works present in literature (in blue) and in this paper (black, green, 
red, and magenta represent results derived from DM, OVISC, CSF, and AGN 
runs, respectively). Thick, and thin lines refer to the total profiles and the DM- 
particle profiles. Finally, the dotted-dashed lines show the regions identified in 
the parameter-space by varying the radial ranges (extreme values in Tabl^2j. 



To conclude this Section, we would like to point that the 
concentration-mass relation derived from the mass profiles of 
the DM particles within the hydro-simulations exhibits normal- 
ization and slope values very similar and close among them and 



to those obtained from the pure N-body simulations. This re- 
sult, joined with the fact that the DM profiles are always well 
represented by an NFW (lower panel of Fig [3]) signifies that the 
best strategy to use in observations would be to compute the 
NFW concentration from the dark matter profile obtained by 
subtracting both the gas and stellar component from the total 
mass profiles. E10 pursed the same idea. They excluded the 
inner 50 kpc to avoid the influence of the stellar component and 
they further subtract the gas mass profile from the mass profile 
derived assuming hydrostatic equilibrium. 

4. OBSERVED C-M 

4. 1 . Synthetic X-ray analysis 

In R12, the X-ray mass profiles have been derived following 
the forward method (Vikhli nTn et al.|[2006| |Meneghetti et al.| 
2010): two fitting formulae have been adopted to fit the sur- 
face brightness and temperature profiles. The analytic best-fits 
have been de-projected, and the mass calculated under the as- 
sumption of spherical symmetry and hydrostatic equilibrium. 
We adopted the mass profiles therein derived to fit the NFW 
formulae using Eq.|2]and Eq.[3] In this setting, the errors on the 
masses are taken directly from the X-ray analysis output (typi- 
cally they span from 10-12% in the most central bins to ~ 5% 
close to R500). 

We remind that the catalogue contains only 20 objects from 
the CSF sample observed at redshift z = 0.25 from three differ- 
ent lines of sight. Therefore, even if the analysis is based on 
60 X-ray images, the sample is not composed of independent 
systems. For this reason and because the mass range covered 
is limited, in the following we will stress more on biases dis- 
covery than on quantifying them thought a precise list of the 
parameters' values. That said, we draw attention to the fact that 
the characteristics of our synthetic catalogue, in number and 
mass range, are similar to the samples of E10 and SA07. 

In Figure [7] we plot three sets of points: 

- c/, Mi are the intrinsic parameters (black asterisks). The 
values are those of the 20 CSF clusters in common with 
the X-ray catalogue and analyzed in the previous Sec- 
tion using the SIM radial range at z = 0.25; 

- Che, Mi w are derived from the 3D mass profiles obtained 
by including in the hydrostatic equilibrium equation the 
true 3D gas density profile and the true 3D mass-weighted 
temperature profile. These values are similarly obtained 
assuming the SIM radial range (black crosses); 

- cx, Mx represent the NFW concentrations and masses 
obtained by fitting the X-ray hydrostatic mass profiles 
with an NFW formula considering the X-RAY radial 
range (red diamonds). However, note that the X-ray 
hydrostatic masses of R12 are an analytic extrapolation 
from 50 kpc to 15% of Rsoq. The latter value indicate 
the boundaries of the central excision applied in R12, 
choice justified by the scope of evaluating the X-ray 
mass bias at large radii. The c-M relation of this set 
of data is shown as a thin red line. The X-ray mor- 



phologically regular clusters as identified in Rasia et al 



(2012a) are shown with the thick red square. The asso 
ciated c-M relation is the thick red line. As found in 
SA07, the NFW formula is generally a good description 
of our X-ray data: only one system presents residual 
above 0.05. This is marked by a magenta square. 
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For direct comparison with observations, we show also the re- 
sults of E10 (blue line). 

Large amount of information can be read from the figure. To 
analyze separately the various levels of complexity involved, 
we will refer also to Figure [8] where we plot two sets of devia- 
tions in both concentrations and masses. 

X — ray — synthetic sample 




[5 



14 M 



FIG. 7. — Concentration and mass values obtained from the 3D intrinsic mass 
profile (cj , M, in black asterisks), the mass profile derived by applying the hy- 
drostatic equilibrium to the 3D gas density and 3D mass-weighted-temperature 
profile (Cke.Mhe in black crosses), the X-ray mass profile of R12 (c x ,M x , red 
diamonds). The c x —M x relation is shown as a red line while ElO's relation is 
reported as a blu line. Thick diamonds and thick red line correspond to mor- 
phologically regular systems. The magenta square indicates the only X-ray 
cluster with a large residual value. 



4.1.1. Deviation caused by the hydrostatic equilibrium 
assumption 

In the top panel of Fig. [8] we show how much concentrations 
and masses vary when we introduce the hydrostatic equilibrium 
assumption in our intrinsic analysis. In this case we define 
A c ( te ,/) = (c he -Ci)/ci and A M(heJ) = (M he -Mi)/M h where i and 
he indicate values shown in Figure |7]by asterisks and crosses, 
respectively. In average, the two masses show a difference of 
25%, being Mj w lower than M, . This value is in agreement with 
the results on the masses presented in R12 (red line in their Fig- 
ure 6), confirming that the NFW fitting procedure is not intro- 
ducing any extra bias. The concentrations show a discrepancy 
of 20% in the other direction without any specific correlation 
with either the mass of the systems, M t , or the mass deviations, 
^M(he,i)- Consequentially, we measure a larger normalization 
for the concentration-mass relation Che-Mi ie when we assume 
the hydrostatic equilibrium but we do not witness any change 
either in slope or in scatter (see also |Lau et al.|2009] ). 

The phenomenon emphasizes that the lack of hydrostatic equi 
librium in each system increases with the radius (Rasia et al 



2004 200 61 |Lauet al.|20 09 Rasi a et al.|2012bl|Khedekaret"aT 
2012 1 Battaglia et al. 2012]l7The profiles of M/, c are closer to the 
true mass profile at the center and diverge more in the external 
region, automatically producing a more concentrated profile. 
The effect, however, is not sensible to the mass of the object, 
as we deduce from the fact that the slope of the concentration- 
mass relation does not vary. This consideration confirms the 
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FIG. 8. — Deviation in concentrations and in masses. Top panel: we compare 
the values obtained from the 3D intrinsic mass profile and the mass derived 
assuming the hydrostatic equilibrium applied to the 3D density profile and 3D 
mass-weighted temperature profile. Bottom panel: these values are compared 
to the X-ray one where we identify also regular and disturbed systems with 
cyan triangles and green diamonds. 



results by Piffaretti & Valdarnini (2008) who, using a more nu- 
merous sample (above 100 objects) and a larger mass range, 
found a negligible dependence on the mass of the hydrostatic- 
mass bias. 

4.1.2. Deviation caused by the X-ray approach 

The outcomes differ when we move to compare the intrinsic 
hydrostatic values (c/,<, and M/ !e ) with the X-ray-derived one (c x 
and M x ). A significant increase of dispersion on both concen- 
trations and masses is present: most clusters at high (low) mass 
have a lower (higher) concentration producing a clear steepen- 
ing of the relation. The deviation A C (x,he) = ( c .v _ Che) I Che = and 
A.M(x.he) = (M x = Mf^/Mhe are shown in the bottom panel of 
Figure [8] 

The majority of the NFW masses from the X-ray data present 
a further reduction, with a median value of -10%. The extra 
drop is caused by the temperature bias affecting this sample: 
the spectroscopic X-ray temperatures are lower than the mass- 
weighted ones because of dishomogeneity in the temperature 
distribution (R12). 

The X-ray concentrations are significantly below the values 
derived from the 3D intrinsic profiles. The median deviation for 
the whole sample is around 35%, and show dependence on the 
X-ray morphology. In Figure [8] we mark with cyan triangles 
and green diamonds the systems belonging to the regular and 
disturbed class of |Rasia et al.| ( |2012a] l (see Appendix of that pa- 
per for the cluster images). The average concentration deviation 
for the sub-sample of disturbed systems reaches 50%. 

The normalization of the c x -M x relation is reduced with re- 
spect to both intrinsics relations. The main reason is the ex- 
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cision (within 15% of R500) applied in R12 to derive the X- 
ray masses. The cut prevents to capture the condensation of 
baryons in the center of the cluster similarly to the process dis- 
cussed above that explains the differences in co when the fit is 
performed in the [0.2-1.27] radial range (Table[2]l. 
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FIG. 9. — Top panel: the concentrations and the masses of the X-ray syn- 
thetic catalogue. Each color is proper of a single cluster observed in the three 
projections. Lower panel: Distribution of the slope of the concentration-mass- 
relation for 30,000 X-ray generalizations. Various percentiles of the distribu- 
tion (5%, 10%, 50%, 90%, 95%) are shown by vertical lines. The thick blue 
line demonstrates that there is less than 3% of probability to obtain a slope 
equal or smaller than the intrinsic value. 



4. 1 .3. Testing degeneracy between parameters 

The steepening of the relation might be a spurious effect orig- 
inated from the fact that the sixty points are not completely 
independent, each single cluster is, instead, represented three 
times. Our concern is instigated by our previous results (Fig. [4] 
and in Fig.[5]l suggesting degeneracy between NFW masses and 
concentrations once they are derived from the similar potential 
wells using slightly different procedural settings (either radial 
range or simulation set or, in this context, separate projections). 

In the top panel of Figure [9] the set of points c x ,M x are plot- 
ted with a choice of color that helps the identification of the 
three projections of the same object. Clearly, this degeneracy is 
present. What needs to be tested is whether the direction of the 



degeneracy is parallel to the relation obtained by using only one 
projection or it is orthogonal and, thus, steepens the relation. 

To test this, we generate 30000 possible combinations of our 
20 objects. In each set, a precise cluster appears only once ob- 
served from a random direction. In this fashion, we are avoid- 
ing the degeneracy caused by having 'same-potential' systems 
in the sample. This approach is more indicated than some re- 
sampling methods such as bootstrap or jackknife since their 
main assumption - independency of the 60 measurements - 
is definitely not satisfied in our case. For each set, we derive 
the concentration-mass relation. The distribution of the 30000 
slopes is shown in the bottom panel of Figure|9] We overplot the 
5 th and 10 th percentiles, the value of the intrinsic slope (black 
asterisks in Figure [7J, and that of the X-ray catalogue (red dia- 
monds in Figure pj. 

The X-ray slope of the c x -M x relation is exactly in the mid- 
dle of the distribution meaning that the most probable slope ob- 
tained by the 20 objects is almost the same to that obtained from 
the enlarged sample of 60 images. From the histogram, we can 
evince also that there is less than 2% probability to reproduce 
the intrinsic slope from randomly orienting each of our clus- 
ters. This test, therefore, reassure that we are not introducing 
an extra bias from our sample. 

4.1.4. Understanding the observational-like c—M relation 

Investigating, further, the reason for the steepening of the re- 
lation we conclude that many aspects, such as projections, dy- 
namical status, dynamical history, hydrostatic-equilibrium mass 
bias, temperature bias, environment, and the dependence of the 
biases on the mass of the systems, are contributing at the same 
time without a clear driving factor. 

Six extreme points are the main responsible to this feature. 
They are identifiable in Figure [7] by having either concentra- 
tion above 10 or mass above 2 x 10 15 h~ l M Q ) (the seventh point 
marked with the magenta square is not accounted to derive the 
relation). All the mentioned six points are related to very mas- 
sive systems. The three with high concentration register a large 
underestimate of the hydrostatic-equilibrium masses and the 
NFW masses, justified by a considerable temperature bias. Since 
this bias increases with the radius, it leads to high concentration 
values. On the other hand, the three points with low concen- 
tration and high NFW mass are three massive and disturbed 
objects located in rich environments. They do not show ma- 
jor substructures but their X-ray emission remains high up to 
the outskirts causing a flattening of the mass profile. Note that 
the X-ray enhancement is present only in one of the three pro- 
jections, therefore there is no intrinsic flattening of the mass 
profile. 

Finally, even excluding these six points, the slope of c x -M x 
is more than lcr larger than the intrinsic slope of c, -M, . The 
explanation is mostly related to two factors noticed above: re- 
stricting the external radial range produces an increase in the 
slope and most massive systems go in both directions of reduc- 
ing the concentration and increasing the mass and vice versa 
(red squares in the bottom panel of Figj4]>. 

5. SELECTION FUNCTION 

Selecting the sample according the X-ray luminosity has been 
recognized as a source of bias in X-ray scaling relations (e.g. 
|Nord et al.|[2008l |Pratt et al.||2"009] l. The wide scatter in the 
luminosity-mass relation induces a concomitant presence of a 
broad range of masses, especially, when small systems are in- 
cluded dPratt et al.|2009j|Stanek et al.|2010l|Mantz et al.|2010[ 
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. The correspondence between the theoretical 
selection - based on mass - and the observational one - based 
on the X-ray emission - is, therefore, compromised. 

At the same time, also the concentration-mass relation presents 
a significant dispersion: galactic or small groups masses can 
have a concentration that vary even by an order of magnitude 
(e.g.|Rudd et al.|2008[|Ludlow et al.|20"l2}|Bahe et al.|2012||Del 
|Boni et al |20TT ~ 



Using the X-ray catalogue we investigate how the two large 
scatters combine together and what is the impact of choosing 
the cluster on the basis of their X-ray luminosity, temperature, 
and SZ-signal instead of their mass. We create various subsam- 
ples selecting systems with: 

- true mass higher than M tlue > lO 15 /! -1 M (subsample Mhigh ) 
and higher than M trlle > 7 x 10 14 /T 1 M Q (M med )\ 

- X-ray luminosity higher than Lx > 14 x 10 44 erg/s (Lhigh) and 
higher than L > 7 x 10 44 erg/s (L med ); 

- SZ-signal higher than SZ > 60 (SZhtgh) and higher than SZ > 

40 (SZ med ); 

- mass-weighted temperature higher than Tmw > 6 keV (Thigh) 
and T MW > 4 keV (T med ). 

The subsamples denoted with the subscript 'high' account for 
either 15 or 18 points, while the 'med' compilations include 35 
-j- 42 points. The mass definition used in this context is based 
on the true 3D mass computed within /?200 while the observ- 
able quantities are measured as projected quantities and inte- 
grated across the whole field of view. The Compton Y param- 
eter called 'SZ' is the integral over the whole field of view of 
adimensional quantity dl/I. 

For each sub-sample we recompute the c-M relation and 
compare the derived slopes. In Figure [10] we color-code the 
different cuts in flux/luminosity as example. Selecting the clus- 
ters by way of their luminosity (subsamples or L med ) con- 
vey the most divergent result with respect to the mass cuts (sub- 
samples Mhigh and M med , respectively). 

The Lhigh clusters show a slope twice as steep as the slope of 
the systems in the Mhigh subsample and 3 times steeper than the 
entire set of objects. Among the smallest clusters, those that 
emit the most have higher concentration. Consequently, a flux- 
limited sample will tend to produce an artificial increase of the 
slope. 

This theoretical finding can be easily translated to an obser- 
vational interpretation: at fixed mass the objects that are more 
X-ray luminous tend to be the morphologically-regular cool- 
core clusters who exhibit a larger concentrations than non cool 
core objects. The X-ray surface-brightness concentration was, 



indeed, the method followed by Santos et al. ( 2008| ) to identify 
cool-core systems at high redshift. The bias toward cool-core 
systems for X-ray cluster sample was also studied in details by 
|Eckert et al.|PDTT) . 

When the objects are selected via their mass-weighted tem- 
perature, we see a good agreement with the c-M relation of the 
mass-selected objects, as expected from the relative small scat- 
ter in the relation Tmw~M. The SZ-selection instead shows a 
slight increase of ~ 15% of the slopes for both SZhi g i, and SZ med 
subsamples with respect to M^gh and M med subsamples. This 
difference is, however, not worrisome being within 1 a. Con- 
cluding our results are suggesting that the SZ selection is more 
similar to the mass (or theoretical) selection, while a compila- 
tion based on the X-ray luminosity substantially contaminate 
the result in the direction of increasing the slope. 

The findings derived from our synthetic catalogue might ex- 



plain the differences among the various observational concentration- 
mass relations. While E10 and SA07, with c-M's slope around 
-0.5, are based on very luminous clusters, the samples of B07 
and |Host & Hansen| ( |201 l| l, with slopes closer to the theoretical 
one, involve also early-type galaxies and small groups that were 
not selected for their X-ray power. The normalization spread 
among these works is also in agreement with our findings about 
the differences caused by the change of the radial range: SA07 
limited the information to very central regions while all the 
other works extended beyond /?2500- 




FIG. 10. — The concentration-mass for the X-ray synthetic catalogue binned 
by the flux emission in units of erg/s. In blue all clusters with X-ray flux (lumi- 
nosity) above 7 X 1CT 12 erg/s (14 X 10 44 erg/s), in cyan those with flux above 
3.5 X 1(T 12 erg/s (7 X 10 44 erg/s), in red all the clusters with an acceptable 
NFW fit. With the same colors, we also plot the best-fit scaling relations as 
solid line and its 1 -sigma variation as shaded region. The magenta squares in 
the bottom of the figure refer to the objects with large NFW residual and not 
used in the c—M relation fitting. 



6. CONCLUSIONS 

Recent claims about the discrepancy between the observed 
and simulated c-M relation were declared real. Here, we tested 
if the difference is instead induced by the two approaches fol- 
lowed to derive the values of concentration and NFW masses. 
Specifically, we checked at the influence of i) the radial range 
adopted to perform the NFW fit of the mass profile, ii) the 
baryonic-physics models included in the simulation, Hi) the proper 
X-ray description, and iv) the impact of the selection function. 
To accomplish the first two tasks we analyzed four sets of 52 
clusters simulated four times varying the description of the ICM 
physics. The latter two aspects were, instead, addressed by 
studying the synthetic catalogue of the 60 Chandra-like X-ray 
images (R12). The four sets of physics considered include dark 
matter (DM), no-radiative runs (OVISC), cooling, star formation 
and feedback by SN (CSF), and cooling, star formation, and 
feedback from AGN (AGN). 

• Radial range. Using the DM set we registered an in- 
crease in both slope and normalization of the c-M re- 
lation when the most external radius considered dur- 
ing the NFW fitting is reduced from 2 x ^200 to 0.4 x 
7?2oo- Extending the fit to inner regions, instead, pro- 
duces an higher normalization accompanied by a shal- 
lower slope. All these differences, even if robust (inter- 
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mediate choices of either internal or external radius re- 
produce the same trend of changes in slopes and normal- 
izations), are mild, Aa ~ 10-20% and Ac ~ 3-5%, 
due to the large scatter of the points. They, however, be- 
come more pronounced once baryons are introduced ( a 
and co increases up to 60% and 15%, respectively). 

NFW fitting and hydro-simulations. The total mass 
profiles of the radiative simulations (CSF and AGN) show 
large residuals in the innermost region . The cause is the 
efficient condensation of baryons in form of stellar com- 
ponent that drastically steepens the total density profile. 
The NR clusters at z = are, instead, well described by 
the NFW formula. For all physics and radial ranges, the 
results related to the DM component are always robustly 
similar. 

Baryon physics. The relations derived from simula- 
tions with baryons have always a larger normalization 
than those from the DM halos. The phenomenon of adi- 
abatic contraction is further aggravated in case of ra- 
diative simulations. Studying the influence of the ra- 
dial range on hydro-simulations we noticed that smaller 
objects are more affected by the change of either in- 
ternal or external radius. The CSF set is characterized 
by a steeper and higher c-M relation (increases around 
30% in both normalization and slope with respect to the 
DM results in the SIM radial range). Both effects are 
strongly reduced when AGNs are included as source of 
feedback. In this case, the slopes are much closer to 
the DM values, with variation of order of few percent. 
Within <10% is also the difference measured from the 
NR set. Comparing the results of our intrinsic analy- 
sis with the observational ones we found that only the 
CSF slope overlaps with the values found in Buot e et aL] 
' 2007[ l. Even in this case, however, it does not reproduce 



the observations by Schmidt & Allen ( 2007); Ettori et al 
( |20T0> . 

X-ray c-M. While the change of the radial range and 
the baryonic physics contribute only minimally to ex- 
plain the gap between the theoretical and observed c—M 
slope, the analysis performed on our X-ray catalogue 
gave more insights. The hydrostatic-equilibrium assump- 
tion introduces a general reduction of the NFW masses 
(25%) and an increase of the concentrations (20%). The 
first result is explained by the HE mass bias. Its ra- 
dial dependence justifies the second finding. The NFW 
masses and concentration derived from the X-ray syn- 
thetic catalogue have a huge dispersion with respect to 
the intrinsic results. Most points show an increase (de- 
crease) of concentration at small (large) masses. The 
slope derived from our X-ray sample is in agreement 
with |Ettori et al.| ( [2010[ ) and is not influence by the con- 
struction of our sample. The relation steepening is caused 
by the concomitant presence of multiple effects such 
as projections, environment, dynamical state, dynami- 
cal history, mass bias, temperature bias, and their de- 
pendences on the radius and on the mass of the system. 
Our morphologically regular objects show a slope 50% 
steeper than the entire sample. This might seem in con- 
trast with the finding of E10 who reported a 20-30% 
shallower relation for their most relaxed systems. How- 
ever, the selection of the two subsamples is based on 



different conditions: in our case regular systems have 
smaller third order power ratio and centroid shift while 
in E10 relaxed objects were defined as having lower 
central entropy level. 

• Selection function By comparing the c-M relations of 
various subsamples defined on the basis of X-ray lumi- 
nosity, mass-weighted temperature, SZ signal, and mass 
we verified the influence of the selection function. The 
clusters selected via the mass-weighted temperature or 
the SZ-signal return similar results than those derived 
from the mass-selected samples. The situation is quite 
different when we choose the systems via their X-ray lu- 
minosity. The c-M slope of the fifteen brightest objects 
is twice as steep as the relation from the fifteen most 
massive clusters. At fixed mass, more luminous clusters 
have systematically higher concentrations. 

Numerical works present in literature are based on many hun- 
dreds of halos. The effects of environment, dynamical status, 
and presence of substructures statistically cancel out when such 
large set of objects is considered. This is not the case when few 
tens of objects are selected as done in observational works. The 
comparison between the two samples, therefore, cannot be as 
straightforward. From our analysis of the intrinsic profile, it 
emerged that smaller objects are more strongly affected by the 
choice of the radial range, and their mass profiles are more eas- 
ily distorted by the mentioned effects. Furthermore, if groups 
or early-type elliptical are chosen for their X-ray power, luck- 
ily, they will not represent the entire population of similar-mass 
objects but they will mostly be characterized by higher concen- 
trations. 

APPENDIX 

In the Table [3] we report the normalizations and slopes of 
the c-M relations derived in the four sets of simulations at 
Z = 0.25 adopting all the radial ranges presented above. The 
DM profiles within the hydro-simulation behave very similarly 
to the profiles from the DM set, therefore, we choose to not list 
them. The significant radii Rsa) and /?200 are recomputed at 
Z = 0.25. Comparing the normalizations of the SIM radial range 
at z = 0.25 with those at z = 0, we notice that they are related to 
each other through: c (z = 0) = c (z = 0.25) x (1 + z)° 4 . 
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Table 3 

Concentration-mass best-fit parameters with 1 a error. The fitting radial ranges are defined in units of R 2 oo (with 

SIM=[0.06-1.27]/?20o) WITH THE ONLY EXCEPTION OF THE X-RAY RADIAL RANGE DEFINED BETWEEN 50 KPC AND 0.6 X ^ 5 00- 
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